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ABSTRACT 

We report on our further analysis of the expanded and revised sample of potential 
BL Lac objects (the 2BL) optically identified from two catalogues of blue-selected 
(UV excess) point sources, the 2dF and 6dF QSO Redshift Surveys (2QZ and 6QZ). 
The 2BL comprises 52 objects with no apparent proper motion, over the magnitude 
range 16.0 < bj < 20.0. Follow-up high signal-to-noise spectra of 36 2BL objects 
and NIR imaging of 18 objects, together with data for 19 2BL objects found in the 
Sloan Digital Sky survey (SDSS), show 17 objects to be stellar, while a further 16 
objects have evidence of weak, broad emission features, although for at least one of 
these the continuum level has clearly varied. Classification of three objects remains 
uncertain, with NIR results indicating a marked reduction in flux as compared to SDSS 
optical magnitudes. Seven objects have neither high signal-to-noise spectra nor NIR 
imaging. Deep radio observations of 26 2BL objects at the VLA resulted in only three 
further radio-detections, however none of the three is classed as a featureless continuum 
object. Seven 2BL objects with a radio detection are confirmed as candidate BL Lac 
objects while one extragalactic (z =0.494) continuum object is undetected at radio 
frequencies. One further radio-undetected object is also a potential BL Lac candidate. 
However it would appear that there is no significant population of radio-quiet BL Lac 
objects. 
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1 INTRODUCTION 

Unification theories of active galactic nuclei AGN) (An- 
tonucci 1993, Urry & Padovani 1995) propose that BL Lac 
objects are the beamed counterparts of FRI/FRII 1 radio 
galaxies. At optical wavelengths BL Lac objects are dom- 
inated by Doppler boosted, synchrotron radiation, the re- 
sultant featureless continuum making BL Lacs all but im- 
possible to find in optical surveys. Until recently, therefore, 
all complete, flux-limited samples of BL Lac objects have 
initially been identified from radio and/or X-ray surveys, 
unification theories seemingly ruling out the possibility of 
these objects being radio- and X-ray-quiet. Seven BL Lac ob- 
jects are detected in the Palomar-Green Survey (PG, Green 
et al. 1986). Three of these were originally classified as DC 

1 double-lobed radio sources classified by Fanaroff & Riley (1974) 



white dwarfs (a continuum white dwarf (WD) with no spec- 
tral absorption features), but were later found by Fleming 
et al. (1993) to be strong X-ray emitters; these authors thus 
consider these sources to be X-ray selected BL Lac objects, 
while the remaining four PG BL Lacs are found in flux lim- 
ited radio samples. 

Identification of the first optically selected sample of BL 
Lac objects from scrutiny of spectra in the 2dF QSO Red- 
shift Survey (2QZ Croom et al. 2001; 2004) was commenced 
in 2000 (Londish et al. 2002, hereinafter Paper I). A similar 
search of colour-selected point sources to optically identify 
BL Lac candidates has now also been conducted using the 
Sloan Digital Sky Survey (Collinge et al. 2005 ). From the 
2dF QSO Redshift Survey we initially selected 56 featureless 
continuum objects as candidate BL Lacs from 6772 spectra 
of objects with 18.25 < 6j < 20.00 and signal-to-noise ra- 
tio (SNR) > 10.0, observed up to August 2001. We have 
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now expanded the survey to include scrutiny of a further 
734 objects satisfying the same selection criteria, observed 
after August 2001 and up to completion of the 2QZ sur- 
vey in April 2002. In addition we have examined spectra of 
bright objects (16.0 < &j < 18.25 and SNR > 10) from the 
6dF QSO Redshift survey (6QZ, Croom et al. 2004); only 
the southern strip of the 6QZ is complete. The survey area 
covered in the 2QZ is 646.8 sq deg and 332.9 sq deg for 6QZ. 
More accurate la proper motion errors in the updated Su- 
perCosmos catalogue 2 led to a further revision of the original 
sample. 

Full details of the selection criteria used (including the 
colour selection of 2QZ and 6QZ objects) are presented in 
paper I and are outlined in part in Section 2 below. In Sec- 
tion 3 we report on results from the high signal-to-noise 
spectroscopy carried out using the Very Large Telescope 3 , 
and in Section 4 on our investigation of the optical and near 
infrared colours of 2BL objects. Follow-up radio observations 
are presented in Section 5. A summary and conclusions fol- 
low in Section 6. 



2 SELECTION OF BL LAC CANDIDATES 

Spectroscopic observations for the 2QZ and 6QZ surveys 
were made using the 2-degree field (2dF) instrument at the 
Anglo- Australian Telescope and the 6-degree field (6dF) in- 
strument at the UK Schmidt Telescope respectively. Data 
reduction and redshift estimation was carried out using the 
2dF/6dF data reduction pipeline (Bailey et al. 2003, Bai- 
ley & Glazebrook 1999). All spectra with signal-to-noise 
ratio (SNR) > 10.0 and 16.0 < bj < 20.00 were visually 
re-examined with the express purpose of identifying objects 
with a featureless spectrum. For consistency we adopted the 
standard operational definition of a BL Lac as an extra- 
galactic object with no emission lines greater than W\ =5A 
and Call H&K break contrast of less than 0.25 (Stocke et al. 
1990). Objects from 2QZ fields with high levels of scattered 
moonlight were excluded from the sample. 

Proper motion studies of the 2BL sample were re- 
peated in June 2003 using the revised SuperCosmos cata- 
logue (Hambly et al. 2001); for the majority of 2dF fields 
a proper motion error of 0.011 arcsec is quoted in each co- 
ordinate, i.e total la error of ~0.015 arcsec. A cutoff of 
2.5a was adopted for inclusion in the 2BL sample, thus we 
are sensitive to transverse proper motions of >45 masyr -1 
for the typical base lines of 11-16 years in the Super Cos- 
mos Sky Survey. The greater accuracy of the June 2003 Su- 
perCosmos catalogue led to the rejection of 15 candidate 
BL Lac objects (see Table 2) included in Paper 1. Large 
la proper motion errors in the original SuperCosmos cata- 
logue in certain fields had resulted in a calculation of < 2.5a" 
proper motion for the 15 objects. One of these 15 objects, 
J140021. 0+001956, was detected at 8.4 GHz (0.336 mjy, see 
§4) and appears to have a nonthermal SED from SDSS pho- 
tometry (Sloan Digital Sky Survey, Stoughton et al. 2002) as 

2 SuperCosmos Sky Survey at www-wfau.roe.ac.uk/sss, main- 
tained by the Institute for Astronomy, Royal Observatory, 
Edinburgh 

3 VLT on Cerro Paranal (Chile) operated by the European 
Southern Observatory, observing run 073.B-0082 
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Figure 1. VLT spectrum of the radio-quiet featureless continuum 
object J215454.3-305654 taken in August 2003. 



well as our own follow-up NIR observations (§4.2). Never- 
theless this object has been removed from the 2BL sample. 
In contrast J113039. 1-004023 which was originally found to 
have > 2.5a proper motion has been re-instated in the sam- 
ple. 

The revised 2BL sample of 52 candidate BL Lac objects 
is listed in Table 1. This updated 2BL sample now includes 
41 objects from the original sample plus the one re-instated 
object, three objects from the 6QZ catalogue and seven ob- 
jects observed with 2dF after August 2001. 

Cross-correlation of this revised sample with the NVSS 
and FIRST radio surveys revealed radio detections for only 
25 per cent of 2BL objects. We therefore conducted our own 
radio observations of the sample, as well as further high 
signal-to-noise spectroscopy and infra-red imaging studies 
in an attempt to rule out contamination of the sample by 
weak-lined AGN and cool, featureless DC white dwarf stars. 



3 HIGH SIGNAL-TO-NOISE SPECTROSCOPY 
3.1 Observations 

It was necessary to obtain high signal-to-noise spectra of all 
2BL objects in order to search for possible weak features not 
distinguishable in the original 2QZ spectra. A spectroscopic 
observing campaign was thus carried out using the FORS1 
instrument at the Very Large Telescope (VLT) with time be- 
ing allocated in service observing mode (program ID 073. B- 
0082) during 2004. To date high SNR spectra of 35 2BL ob- 
jects have been obtained at the VLT. J215454. 3— 305654 was 
first observed as a test case in August 2003 (during observing 
run 71.A-0174) and confirmed as an extragalactic feature- 
less continuum object at z = 0.494 from weak absorption 
features (CaH&K and G-band) and [OIII] emission (Figure 
1; see Londish et al. 2004 for full details). Spectroscopy of 
the remaining 34 objects was obtained during 2004 (April- 
August), using the long-slit spectrograph (larsec slit) and 
300V grism with central wavelength 5900 A. Exposure times 
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Table 1. 2BL sample of 49 objects from the 2QZ survey and 3 objects from the 6QZ survey 
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* These objects have been added to the original sample presented in Paper I;** signifies an object from the 6QZ survey. 

a 1.4 GHz fluxes are from NVSS (NRAO/VLA Sky Survey, Condon ct al. 1998) and/or FIRST (Faint Images of the Radio Sky at 20cm, 

White et al. 1997) and upper limits from Australia Telescope Compact Array measurements (see §5); 8.4 GHz fluxes and upper limits 

are from our own VLA observations (§5); 843 MHz flux from the SUMSS Catalogue (Bock ct al. 1999). 

6 AGN I: broad/weak line AGN 

c BLL signifies that the object is classified in the literature as a BL Lac (Vcron-Cetty & Veron 2001; Voges et al. 1999) 
fOptical photometry from SDSS, NIR photometry from Kitt Peak, Nordic Optical Telescope and Calar Alto 

ft Overall classification from our VLT observations (§3) and/or SDSS (Sloan Digital Sky Survey, Stoughton et al. 2002) spectra (Fig 
Al), and NIR imaging (§4). 

Redshifts are as quoted in SDSS and/or measured from VLT spectra - see Figs 3-5 
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Table 2. The 15 featureless continuum objects removed from the 2BL in June 2003 on account of their > 2.5c proper motion. Objects 
marked f have SDSS photometry (see §4) that indicates the object is indeed most likely stellar, whereas for the two objects marked ft 
no blackbody curve could be fit to the SDSS optical magnitudes 
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varied from 2 x 10 to 2 x 15 mins, depending on the bright- 
ness of the source. The observations were made in generally 
poor seeing conditions (0.7"-2.2") and six objects were ob- 
served during bright time. Flux standards were not observed 
on all nights. Flux calibration was nevertheless performed 
using standards taken on another night in order to correct 
for instrument response. This, together with slit losses, has 
resulted in the flux calibration of our 2BL objects being only 
approximate, with spectra showing a drop off in flux below 
~4200 A; contamination from second-order overlap redward 
of 7500 A is also evident in the spectra of blue objects. 

Data reduction was carried out using standard IRAF 4 
routines. These high signal-to-noise spectra (Figs. 3—5) show 
features not clearly defined in the original 2QZ spectra. 



fit and to discriminate between the two models, so the exact 
value of the error Oi is not crucial. However, if the reduced- 
X 2 value (defined as xt = X 2 l v where v is the number of 
degrees of freedom, in this case two less than the number 
of data points) is greater than 5, the fits are rejected. In 
four cases where emission and absorption features dominate 
the spectrum neither a power law nor a blackbody function 
could be fit. In all other cases the fit clearly indicated either 
a blackbody curve or positive powerlaw slope (i.e increasing 
flux for increasing wavelength). The actual value of T or a 
returned by the fitting routine will differ depending on the 
choice of regions excised from the fitting procedure. For ther- 
mal sources the estimated la error on the best fit blackbody 
temperature is 100-250 K. 



3.2 Curve fitting of VLT spectra 

In order to classify 2BL objects with VLT spectra as either 
thermal or nonthermal a curve fitting procedure was per- 
formed using a least-squares fit to the data. A list of results 
is shown in Table 3. For each spectrum, two models are fit, 
each with two free parameters: a blackbody function (tem- 
perature and normalisation) and a power law function (slope 
and normalisation). Fitting was performed over the wave- 
length range ~4000-8000A such that instrumental effects do 
not contribute to the fit. In addition those parts of the spec- 
trum strongly affected by features such as atmospheric lines 
or emission and absorption features are blanked out interac- 
tively. In this way the models are fit to the estimated con- 
tinuum via a least-squares (or x 2 -niinimisation) technique. 
A fixed error on each data point of Oi = 2/iJy was assumed. 
The model with the lower value of \ 2 = — m 2 /°f 

(where /; are the measured fluxes, and nii are the model 
fluxes) is chosen as the best fitting model. Note that the 
value of x 2 is used only to ascertain a relative goodness of 

4 Image Reduction and Analysis Facility supported by the Na- 
tional Optical Astronomy Observatories (NOAO) in Tucson, Ari- 
zona. The NOAO is operated by the Association of Universi- 
ties for Research in Astronomy (AURA), Inc. under cooperative 
agreement with the National Science Foundation. 



3.3 Results 

Spectra of 14 objects appear to be well fit by a thermal 
blackbody curve, although one of these objects - J224559.1- 
312223 - has a weak emission feature that could be Mgll 
at z — 1.324. Of the remaining 21 objects seven appear 
to be featureless continuum objects (only one of which, 
J215454. 3— 305654, is radio-quiet), and 13 show weak, often 
very broad emission features as well as absorption systems 
in their spectra. One of these 13 objects, J014310. 1—320056, 
is classified in the literature as a BL Lac, and does indeed 
appear to have a featureless continuum in the 2QZ spectrum 
initially used to select the object as a candidate BL Lac (Fig 
2, left). However in a second 2QZ spectrum taken 9 months 
later (Fig 2, right), and in our VLT spectrum, the object is 
clearly a broad line AGN (z = 0.375). The spectrum of the 
remaining object (J113413. 4+001041) is too noisy to distin- 
guish the presence or absence of lines; this object is classed 
as a QSO at z = 1.4856 from its SDSS spectrum (Fig Al). 

In addition to J113413.4+001041 another eight objects 
now also have spectra available in the Sloan Digital Sky 
Survey Data Release 4 (SDSS, Stoughton et al. 2002). Five 
of these flux calibrated SDSS spectra (Fig Al) confirm the 
identity of two objects as candidate BL Lacs, two as QSOs 
and one as a WD (white dwarf) star; SDSS spectra of a 
further three 2BL objects without VLT spectra identify the 
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Figure 2. 2QZ spectra of J014310. 1-320056 taken on 12/01/00 and 17/09/01. Emission features in the spectrum on the right put the 
object at z = 0.3751 



objects as QSOs. Classification of all 2BL objects is listed 
in Table 1. 

VLT spectroscopy indicates ~40 per cent contamination 
of the sample by stellar objects with a similar percentage of 
contamination by AGN with weak or broad features. In the 
case of the latter it is not possible to determine from the 
original 2QZ/6QZ spectra (which are not flux calibrated) 
whether objects were misclassified or whether continuum 
levels have subsequently changed. A comparison of origi- 
nal 2QZ/6QZ spectra (corrected for 2dF/6dF instrument 
response) with the corresponding VLT spectrum is shown 
in Appendix C. 

Note that one radio-loud 2BL source without a VLT 
spectrum (J225453.2-272509) was observed in 2001 at the 
Siding Spring Observatory, Australia (Fig 6 - see Paper 1 
for details); it is classified as a continuum source at a red- 
shift of z =0.333 identified from weak Call H&K absorp- 
tion features and G-band in the spectrum. Paper 1 also in- 
cludes high signal to noise Keck spectra of the two radio- loud 
objects J105534.3-012617 and J121834.8-011955 (observed 
by Brotherton et al. 1998), showing these objects to be fea- 
tureless continuum sources. 



J225453. 2-272509 
1.5xl0->« i i i i | i i i i | i i i i | i i i 




4000 4500 5000 5500 6000 



MA) 



Figure 6. J225453.2-272509, a 52mjy continuum object at 
z =0.333 



4 BROADBAND OPTICAL AND NEAR 
INFRARED STUDIES 

4.1 SDSS Magnitudes 

Prior to obtaining the above VLT spectra we also checked 
the Sloan Digital Sky Survey for photometry of candidate 
BL Lac objects. SDSS Data Release 4 now includes pho- 
tometric magnitudes 5 for 19 2BL objects in the equatorial 



strip of the 2QZ (the only exception being 2QZJ105534.3- 
012617). 

Curve fitting was also carried out on the 19 2BL ob- 
jects with SDSS optical magnitudes; good fits to a black- 
body function were found for five objects, while a further 
three were better fit by a blackbody than a powerlaw slope. 
Of the remaining 11 objects three were found to be well fit 
by a powerlaw slope (J113900.5-020140, J121834.8-011955 
and J142526. 2— 011826), one was better fit by a powerlaw 
than a blackbody function (J103607. 4+015658) and seven 



5 SDSS uses photometric magnitudes designated v! , g' , r', i' and 
z', with central wavelengths of 3500 A, 4800 A, 6250 A, 7700 



A and 9100 A respectively, and AB magnitude system with a 
common zero point for all filters (Fukugita et al. 1996) 
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Figure 3. VLT spectra of 2BL objects, fitted with either a blackbody curve or power-law slope. In three cases no slope could be fitted. 
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Figure 4. VLT spectra of 2BL objects, fitted with either a blackbody curve or power-law slope. 
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Figure 5. VLT spectra of 2BL objects, fitted with either a blackbody curve or power-law slope. In one case where absorption and emission 
features dominate the continuum no slope has been fitted 
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Table 3. Curve fitting results for 2BL objects with VLT spectra and SDSS optical photometry 



object 






UL^lllLJ I (_X 

Tynivi \/ 1 * 1 QT^^r^TT 1 ^ 

llt_Jlll V 1_J _1 o L/L^t-* LI a 


UClllp / ( t 
11 Will ulVkJu L/11W LrWlllC \iL y 


J002522.8-284034 


0.00 


Odd WD with split lines 


25,500 K 




J002746.6-293308 


0.00 


WD with He and Ho 


11,600 K 




J00275L 5-293506 


0.00 


no obvious features 


10,400 K 




J003058.2-275629 


1.606 


Broad Mgll, CIII] 


0.37 




J004f06.8-291ff4 


0.00 


DB with Hel and Ha 


10,100 K 




J004950.6-284907 


0.00 


Weak Balmer lines 


6,900 K 




J0143f0.f-320056 


0.375 


Strong Balmer lines, Mgll (Si. 4 =76mjy) 


1.1/1.85 




J022044.7-3037f3 


1.439 


Weak, broad Mgll, CIII], abs at z=1.202 


- 




J023536.7-293843 


- 


BL Lac obj 


0.70 




J0304f6.3-2832f7 


- 


BL Lac obj 


1.14 




J031056.9-305901 


1.911 


Broad Mgll, weak CIV, CIII] 


- 




Jf03607.4+0f5658 


1.875 


Broad Mgll, CIV, CIII] 


0.71 


0.43 


Jf05355.f-005538 


0.00 


Weak Balmer lines 


8,700 K 


9,500 K 


J105534.3-0f2617 


- 


BL Lac obj 


1.1 




Jfl0644.5+0007f7 


— 


- 


- 


27,000 K 


Jlf3413.4+00f04f 




Noisy spectrum 


- 




Jlf3900.5-020f40 


- 


- 


— 


1.57 


Jfl4137.f-002730 


- 


- 


- 


0.67 


J114221.4-014812 


- 


- 


- 


10,800 K 


Jfl4554.8+00f023 


- 


- 


- 


6,800 K 


Jf200f5.3+000552 


1.650 


Broad Mgll, weak CIV, CIII] 


2.1 




J121834.8-0fl955 


- 


strong continuum object (Si. 4 =244mjy) 


1.74 


1.44 


Jf23437.6-0f2953 


1.716 


Broad CIV, CIII], Mgll 


1.25 




J125435.7-0fl822 


0.00 


no obvious features 


6,800 K 


7,200 K 


J125501. 2+015513 


- 


- 


- 


10,400 K 


J140207.7-0f3033 


0.00 


DA 


14,500 K 


13,400 K 


Jf4f040.2-023020 


0.00 


no obvious features 


9,100 K 


10,500 K 


Jf42526.2-0fl826 


- 


continuum object (Si. 4 =10 mjy) 


1.6 


1.12 


J2f5454.3-305654 


0.494 


weak [OIII], Ca H & K, G-band 


0.36 




J220222. 7— 285306 


1.231 


ttt 11 1 -K r TT 1 /~1TTT1 1 T~l TT 

Weak broad Mgll and CIII] plus Fell 


1.15 




J220850.0-3028f7 


1.711 


Weak broad Mgll, CIII] and CIV plus Fell 


1.03 




J221f05.2-284933 


1.840 


Weak broad Mgll, CIII] and CIV (S1.4 =81mjy) 


0.7 




J221450.f-293225 


1.626 


Weak broad Mgll, CIII] and CIV plus Fell 


1.55 




J223233.5-272859 


0.00 


DB 


13,500 K 




J224559.f-3f2223 


?? 


star or Mgll at z=1.324 


7,100 K 




J230306.0-312737 


2.459 


Weak broad CIII], CIV and Ly-alpha 


1.28 




J230443.6-311107 


0.00 


WD, DO, Weak Hell at 4686 A 


10,000 K 




J231749.0-285350 


0.00 


WD, Ho + Hel lines 


10,400 K 




J232531. 3-313136 


0.00 


WD, weak Ho, H/3 


9,500 K 




J233508.1-283035 


1.386 


Very weak broad Mgll and CIII] (Si. 4 =6mjy) 






J234414.7-312304 




continuum object (Si. 4 =5.4mjy) 


1.86 





WD classifications: DA denotes a WD with a hydrogen-rich atmosphere or outer layer, indicated by strong Balmer hydrogen spectral 
lines; DB, a neutral helium-rich atmosphere indicated by neutral helium spectral lines, (He I lines); DO, an ionized helium-rich 
atmosphere indicated by ionized helium spectral lines, (He II lines); DC, no strong spectral lines. 

Objects have been fitted with cither a blackbody curve or a power-law slope characterised by S oc v~ a . The best fit value of T or a is 
shown. See §3.2 and §4.1 for details. 



produced a poor fit in both cases. Fits to optical and NIR 
photometric magnitudes are listed in table 3 and shown in 
Appendix B; the unusual SED of J113039. 1-004023 is a re- 
sult of broad absorption and emission features, evident in 
the SDSS spectrum of this object. 



4.2 Infrared imaging 



searched the Two-Micron All-Sky Survey (2MASS) 6 for NIR 
magnitudes of our 2BL sources. We found only one detec- 
tion of the bright 6QZ object J022044. 7-303713 which has 
a J — K colour of 0.67, indicative of a nonthermal source. 
As 2MASS limiting magnitudes for point sources are 15.8, 
15.1 and 14.3 for the J- H- and ^ s -bands (1.25, 1.65 and 
2.17 microns respectively), no 2QZ objects are found in the 
2MASS catalogue, although J105534.3-012617 has a 3cr de- 
tection. We therefore carried out our own infrared imaging 
as a further means of distinguishing between cool, feature- 



Given that the SEDs of both thermal and nonthermal ob- 
jects can be similar at optical wavelengths (see Fig. 7) wc 



6 http: / / www.ipac.caltech.edu/2mass / 
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LHS 2333 (White Dwarf) 
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Figure 7. The energy distributions of 2 DC WDs (T e ff ~ 
7000K) from Bergeron ct al. (2001) (solid symbols) compared to 
that of a luminous (V=16.4) BL Lac (asterisks) (data from Tanzi 
et al. 1988) and an / <x u~ 0,5 energy distribution typical of an 
AGN/QSO. For the 2 WDs the J - K colours are +0.19 (G24-9) 
and +0.12 (LHS 2333). Note the 3-4 leftmost points arc insuf- 
ficient to differentiate between the thermal, blackbody curve of 
WDs and the nonthermal continuum of an AGN; thus in some 
cases near-IR data points are essential to determine the fit. 



less DC white dwarfs and genuine candidate BL Lac objects. 

NIR imaging was first carried out in April/May 2002 
using the infrared camera, SQIID, on the 2.1m telescope 
at Kitt Peak, Arizona. This instrument produces simultane- 
ous images in the J, H , K s , and narrowband L passbands. 
With 440 x 460 pixels and 0.69 arcsec per pixel, the field of 
view is 304 x 317 arcsec (~ 5 x 5 arcmin). Candidate BL 
Lac objects were observed in J, H and K (corresponding 
to an effective central wavelength of 1.267, 1.672 and 2.225 
microns respectively). 

Sources were observed using 20 co-added 15s exposures 
and a 5-image dither pattern. Spectroscopic standards were 
chosen from Persson et al. (1998). Objects and standards 
observed on the first night, which was not photometric, were 
briefly re-observed at the beginning of the third night. Data 
reduction was performed using standard IRAF routines. 

Of the six 2BL objects observed, four were subsequently 
removed from the sample on account of their > 2.5<r proper 
motion. All four were found to have thermal NIR mag- 
nitudes, consistent with their optical SDSS magnitudes. 
Optical/IR SEDs obtained for the two remaining objects, 
J113900.5-020140 and J123437.6-012953 are shown in Ta- 
ble 4. Although the SED of these objects cannot be well fit 
by a single powerlaw slope, their J — K values are consistent 
with those found for QSOs by Barkhouse and Hall (2001) 
using 2MASS sources. Similar NIR magnitudes were found 
for these objects in later observations at the NOT and Calar 
Alto (sec Fig. Bl). 

Infrared observations of 17 2BL objects were also car- 
ried out during April 2003, January 2004 and March 2004 



on the 2.56m Nordic Optical Telescope , La Palma, Canary 
Islands, using the wide field camera (NOTCam) providing 
a field of view of 3.9' x 3.9' (0.23" /pixel). J, H and K fil- 
ters at 1.235, 1.662 and 2.159 microns were used. Results 
are presented in Table 4. Data reduction was carried out 
with IRAF's APPHOT package using standard NIR tech- 
niques (sky-subtraction using a sky-frame evaluated from 
a 12-frame set of dithered images, followed by Hatfield divi- 
sion using twilight flats, and finally shifting and combining). 
With the exception of one night (Jan 2004) NOT data was 
taken in non-photometric conditions, hence data quality is 
not very high. The target magnitudes were estimated using 
differential photometry against 2MASS field stars. Typically 
more than two 2MASS stars were available. As a check, data 
from the 1st of January 2004 was calibrated using standard 
stars from Hunt et al. 1998 and no significant differences 
were found between the 2MASS and Hunt et al. calibrations. 

Three 2BL objects now removed from the sample on 
account of their high proper motion were also observed; two 
are confirmed to be thermal however from both SDSS and 
NOT magnitudes the source J140021.0+001956 would ap- 
pear to be a nonthermal object. A 3-ct radio detection of 
0.113 mJy was also recorded for this object with the VLA 
(§3, Table 5). 

Further NIR images in J and/or K' (1.250 and 2.160 
microns respectively) of 17 2BL objects were taken at the 
Calar Alto 2.2m telescope 8 during the period May, 19-27 
2005. The NIR camera MAGIC in high-resolution mode 
was used, which provides a field of view of ~ 2.7' x 2.7' 
(0.636" /pixel) by means of a 256 x 256 pixel HgCdTe array. 
For the observations, we adopted a 3 x 3 dither pattern, with 
the individual raster points separated by 10" . Depending on 
the filter, at each of the individual raster points 15 images 
2sec each in J or 30 images lsec each in K' were recorded, 
summed and saved. The total integration time for one dither 
pattern was thus 450 sec. Since our aim was to detect the 
source and measure its (unknown) flux, the integration times 
in total per filter and object were adopted adaptively. If the 
quick online reduction at the telescope shows the object to 
be detected after one dither pattern, the observations of this 
particular object were either repeated using the other filter 
or stopped. Otherwise the observations in a particular fil- 
ter were repeated until the object was detected on the fi- 
nal summed image. If the object was not detected after 4 
dither patterns (~ 30min), the observations of this source 
were terminated as well. Observing conditions varied greatly 
during our 8 night run (4 full and 4 half nights), and obser- 
vations taken during poor seeing conditions (>2" FWHM) 
were repeated. Only two of the nights were photometric dur- 
ing which we regularly observed standard stars from the 
lists of Elias et al. (1982) and Hawarden et al. (2001) to set 
the zero points. The data reduction was standard for NIR- 
observations (subtracting a scaled sky derived from adja- 

7 observations made with the Nordic Optical Telescope, operated 
on the island of La Palma jointly by Denmark, Finland, Iceland, 
Norway, and Sweden, in the Spanish Obscrvatorio del Roque de 
los Muchachos of the Instituto dc Astrofisica de Canarias 

8 observations collected at the Centro Astronomico Hispano 
Alcman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck Institut fur Astronomie and the Instituto de Astrofisica 
dc Andaluci'a (CSIC) 
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Object 


J 


Kitt Peak 
H K s 




Classification 


J113900.5-020140 
J123437.6-012953 


17.78 ±0.06 
18.58±0.09 


17.45 ±0.07 16.94 ±0.07 
17.99 ±0.10 17.88 ±0.11 




?? 
?? 


Object 


J 


Nordic Opt Telescope 

H K 


Calar Alto 
J K' 


Classification 



J105355. 


1-005538 


19.111 ±0.092 




> 


19.0 












stellar 


J105534. 


3-012617 


16.884 ±0.087 


16.237 ±0.067 


15.540 ±0.130 


16.97 ±0.22 


15. 


70 


±0. 


16 


nonthermal 






16.739 ±0.094 


16.069 ±0.055 


15.280 ±0.118 














J110644, 


5±000717 


> 20.0 


> 19.2 


> 17.7 






> 


17 




?? 


J113039. 


1-004023 


18.451±0.081 


18.340 ±0.124 


17.734 


±0.111 




18. 


.27 


±0. 


31 


?? 


T1 1 9/11 Q 

J llo41o. 


4-)-UUlU41 












i *7 k*7 _l_n 1 1 


16 


.80 


±0. 


17 


77 


J113900 


5-020140 


17.848 


±0.082 


17.771 ±0.119 


17.241 


±0.119 


18.24 ±0.20 


17. 


20 


±0. 


21 


variable? 






17.617 


±0.099 


17.191 ±0.080 


16.955 


±0.150 














J114137. 


1-002730 


18.199 


±0.055 


17.837±0.085 


17.768 


±0.174 


17.66 ±0.15 


16. 


.65 


±0. 


15 


variable? 


J114221. 


4-014812 


17.739 


±0.122 


17.353 ±0.086 


16.884 


±0.118 


17.70 ±0.26 


17. 


26 


±0. 


24 


stellar 


J114327. 


3-005050 


18.294 


±0.098 


17.480 ±0.131 






18.48 ±0.34 


17. 


49 


±0. 


22 


nonthermal 






18.112 


±0.128 


17.511 ±0.131 


17.137 


±0.116 














J114554. 


8±001023 


19.216 


±0.077 




19.683 


±0.270 






> 


18 




stellar 


J115909. 


6-024534 


18.193 


±0.079 




17.467 


±0.106 


18.76 ±0.46 


17. 


32 


±0. 


20 


variable? 






18.305 


±0.083 




17.538 


±0.106 














J120015. 


3±000552 


17.973 


±0.058 


17.268 ±0.093 


17.145 


±0.054 












?? 


J121834. 


8-011955 


15.496 


±0.136 


14.543 ±0.139 


13.574 


±0.049 












nonthermal 


J123437. 


6-012953 


18.521 


±0.098 




17.470 


±0.114 




17. 


.58 


±0. 


25 


?? 










17.990 ±0.067 


17.517 


±0.103 














J125435. 


7-011822 


18.733 


±0.087 


18.557 ±0.198 


18.167 


±0.168 




17. 


.63 


±0. 


54 


stellar 


J140207. 


7-013033 


> 19.3 ±0.5 




> 18.2 ±0.5 






> 


17.3 






J141040. 


2-023020 


19.239 


±-0.232 


















stellar 


J142526. 


2-011826 


17.656 ±0.164 




15.972 ±0.127 


18.13 ±0.25 


16.15 ±0.08 


nonthermal 






17.793 ±0.093 





















Objects with high proper motion now removed from the 2BL sample 



J102615. 3-000629 
J130009. 9-022600 
J140021.0+001956 



18.988 ±0.132 
18.129 ±0.176 
18.494 ±0.082 



19.097 ±0.210 
18.427 ±0.211 
17.508 ±0.055 



stellar 
stellar 
nonthermal 



Table 4. J, H and K/K' magnitudes of objects observed at Kitt Peak, the Nordic Optical Telescope and Calar Alto 



cent observations of the same field, flatfielding by means of 
a dome flat taken as the difference of exposures of the illu- 
minated and non- illuminated dome). Bad pixels were iden- 
tified from the flatfield exposures and subsequently flagged 
on the individual images. The reduced images were then 
aligned to integer pixels and averaged by a clipping algo- 
rithm to get rid of the bad pixels. Photometry was done 
using SExtractor (Bertin & Arnouts 1996)), where we used 
the photometry provided by Mag_Best. As only a few nights 
were photometric calibration for all sources was performed 
via 2MASS stars in the same frame. The zero points derived 
via 2MASS observations and those derived from our own 
standard star measurements agree very well. The resulting 
magnitudes and errors for our 2BL objects are given in Ta- 
ble 4. As can be seen, the errors are partly very large. This 
is due to the fact that for the majority of our sources only 
1-2 2MASS stars were present in the same frame, and these 
stars typically had large measurement errors. Thus our er- 
rors are dominated by systematic errors. The internal errors 
are much smaller - of the order of 0.02-0.05mag. 



5 RADIO PROPERTIES 

Only nine of the 2QZ BL Lac candidates have catalogued 
FIRST and/or NVSS detections. However, because of their 
faint optical flux densities (mj, > 18.25), it is possible that 
some of these objects are moderately radio-loud BL Lacs at 
higher redshift that simply fall below the flux limits of these 
surveys. Prior surveys have found that BL Lacs, and other 
radio-loud AGN, tend to lie above a ro > 0.2 whereas radio- 
quiet AGN lie below a ro < 0.2 (Rector et al. 2000). The 
NVSS and FIRST survey detection limits are insufficient to 
determine whether or not the 2QZ BL Lac candidates are 
radio-loud objects at higher redshift or if they are indeed 
radio-quiet BL Lacs that occupy a distinct region of the 
olox, a ro plane. 

Twenty-seven objects from the 2QZ sample that lacked 
detection in the FIRST and NVSS surveys were observed 
with the NRAO 9 Very Large Array. The D configuration 



9 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated by Associated Universities, 
Inc. under a cooperative agreement. 
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Object 


Sg.4 (mjy) 


rms 


J114137. 1-002730 


3.580 


0.052 


J115909. 6-024534 


0.296 


0.037 


T1 40021 0+001055+ 


0.336 


0.034 


J031056. 9-305901 


0.129 


0.035 


T140207 7-0130331 + 


0.113 


0.038 


J004950.6-284907 


<0.13 


0.042 


J023405.5-301519 


<0.10 


0.032 


J103607.4+015658 


<0.11 


0.038 


J105355. 1-005538 


<0.14 


0.048 


J110644.5+000717 


<0.12 


0.038 


J113413.4+000411 


<0.21 


0.042 


J113900.5-020140 


<0.10 


0.033 


J114221.4-014812 


<0.11 


0.035 


J114327.3-005050 


<0.11 


0.037 


J120015. 3+000552 


<0.14 


0.045 


J120558. 1-004215* 


<0.12 


0.039 


T100QQO n n 1 t:ci(i* 


<0.11 


0.036 


J123437.6-012953 


<0.08 


0.026 


J125435. 7-011822 


<0.11 


0.035 


J130009.8-022600* 


<0.22 


0.073 


J132811. 5+000227* 


<0.09 


0.031 


J131635. 1-002810* 


<0.11 


0.035 


J223233.5-272859 


<0.14 


0.047 


J224559. 1-312223 


<0.13 


0.044 


J230443. 6-311107 


<0.13 


0.043 


J231749. 0-285350 


<0.12 


0.039 


J232531. 3-313136 


<0.14 


0.048 



Table 5. 8.4 GHz flux densities (or upper limits) of 23 2BL 
sources, f This object is assigned a proper motion of 3.87(7 in the 
updated SuperCosmos catalogue and has been removed from the 
sample. However IR observations conducted at the NOT identify 
this object as a nonthermal source, consistent with SDSS opti- 
cal photometry. * These objects have now been removed from the 
sample on account of their > 2.5(7 proper motion, ff This object 
is identified from the VLT spectrum as an unambiguous WD. 



was chosen to maximize phase stability. The targets were 
observed at a frequency of 8.4 GHz on 20 and 21 April 2003, 
with all 27 antennas operating. Each object was observed for 
14 to 16 minutes. The data were reduced using the NRAO 
AIPS software package. Flux density calibrations were boot- 
strapped from 3C 48 on April 20th and 3C 286 on April 21st 
respectively. All objects within the primary beam were im- 
aged and cleaned to remove residual sidelobes from nearby 
bright radio sources. The noise levels in the resultant maps 
were a s ^ 30-50 /iJy beam -1 . 

Of the 27 objects observed, only four detections were 
made coincident with the optical positions of the BL Lac 
candidates. All other sources were undetected to within flux 
densities of S pea k > 3a above the background noise level. 
The VLA observations are summarized in Table 5. Column 
1 is the 2QZ object name. Columns 2 and 3 give the peak 
detection and noise level for each object in milli-Janskys. 
For non-detections, the peak value is calculated to be a 3a 
upper limit. 

Two objects were also observed at 1.4 GHz using Di- 
rector's time at the Australia Telescope Compact Array 
(ATCA), Narrabri NSW, in November 2003 and May 2004. 
No detection for either J215454.3-305654 or J003058.2- 



275629 was found, placing a 3a upper limit on the flux den- 
sity of these sources of 0.21mJy and 0.15mJy respectively. 

Radio detections have thus been found for 15 of the 
candidate BL Lacs. 



6 SUMMARY AND CONCLUSIONS 

From VLT spectra of 35 2BL objects we identify 11 WDs 
with weak or absent absorption features over a temperature 
range of 6800 - 13500 K, one unambiguous DA at 14500 K, 
plus a further unusual WD with split absorption features 
at a temperature of ~ 25000K. This represents contamina- 
tion of the 2BL sample by stellar objects of ~40 per cent. 
NIR observations at the NOT and Calar Alto identify a 
further two stellar objects, as well as three objects clas- 
sified as stellar from the VLT spectra. Two other objects 
(J110644.5+000717 and J125501. 2+015513) have SDSS op- 
tical magnitudes that are well fit by a thermal blackbody 
curve. J110644.5+000717 was also observed at the NOT and 
Calar Alto but was too faint for a detection. NIR upper lim- 
its are consistent with the SDSS magnitudes that identify it 
as a thermal object. Note however that Nesci et al. (2005) 
report this object as showing evidence of variability. 

Twelve objects with VLT spectra (four of which arc 
radio-detected) are identified as weak-lined AGN; eleven of 
these are in the redshift range z = 1.23-1.91 where the Mgll 
line has fallen in the noisy red end of the 2QZ spectra. A 
further X-ray- and radio-loud object (J014310. 1-320056) at 
z =0.375, has clearly varied since the first 2QZ spectrum was 
taken. Nine 2BL objects have SDSS spectra (Fig Al); three 
are of objects without VLT spectra, and all three appear to 
be AGN. 

Ten 2BL sources remain without a clear classification. 
For three of these objects with SDSS and NIR photome- 
try it would appear that either the source has varied con- 
siderably between the NIR observations and those carried 
out for SDSS or these are a new class of object with an 
unusual kink in the spectrum between the optical and in- 
frared. Note that within the errors the NIR observations 
carried out at the NOT in 2003 and 2004 and at Calar Alto 
in 2005 are generally in agreement. One of these three ob- 
jects, J114137.1-002730, is quoted by Nesci et al. (2005) as 
having high proper motion, but would appear from our ob- 
servations to be a 3.58mJy radio source. For the other seven 
objects only the 2dF/6dF spectra are available. 

A total of eight 2BL objects are classed as continuum 
objects. Seven are radio-loud and three of these are detected 
at both X-ray and radio frequencies. NIR observations con- 
firm the nonthermal nature of three of these eight continuum 
objects. J215454. 3— 305654 is undetected at radio frequen- 
cies. The VLT spectrum confirms that it is extragalactic, at 
a redshift of 0.494. As the spectrum does show weak, nar- 
row [OIII] emission this object could belong to the recently 
discovered population of low luminosity AGN with no broad 
line region (Nicastro et al. 2003; Elitzur & Shlosman 2006). 
This object has been included in a Chandra proposal to 
obtain X-ray data for a sample of radio- weak/radio-quiet 
featureless continuum objects (PI J. Gelbord, MIT). A fur- 
ther object, J224559. 1-312223, is potentially a radio-quiet 
BL Lac; the VLT spectrum shows a weak feature that could 
be Mgll at z =1.324, however the shape of the spectrum 
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Table 6. Classification of 2BL sources from VLT spectra, NIR 
imaging and SDSS spectra and photometry 



number 



classification 



17 
16 
10 
7 
2 



stellar 
AGN 

?? 

continuum (radio-detected) 
continuum (no radio detection) 



is also well fit by a blackbody curve of T = 7100K, and the 
feature could then be Ha at zero redshift. A colour-colour 
diagram of these nine 2BL candidate BL Lacs is shown in Fig 
8 (Top); the two radio quiet candidates lie well outside the 
region of colourspace populated by hitherto discovered BL 
Lacs. These two objects are now the only remaining radio- 
quiet candidate BL Lacs of 47 objects included in Paper 1 
(see Fig. 8 (Bottom)). 

A similar search for optically selected candidate BL Lac 
objects (or QFOs, quasi-featureless objects) was carried out 
using the SDSS survey (Collinge et al. 2005). This search 
made use of colour-selection to separate extragalactic ob- 
jects from galactic QFOs - mostly DC white dwarfs - where 
no reliable proper motion measurements were available. Of 
240 probable QFOs identified from the 2860 sq deg of the 
SDSS spectroscopic survey, 27 objects are currently unde- 
tected at the 5cr level of the FIRST and NVSS surveys, how- 
ever Collinge et al. (2005) indicate that more sensitive radio 
observations are needed before deeming these objects to be 
radio-quiet. In comparison the 2BL sample, from a survey 
area of ^650 sq deg and magnitude range 18.25 < 6j < 20.0 
(i.e. not including the 6QZ survey as all three candidates are 
now classified as AGN), identifies seven BL Lac objects plus 
two radio-quiet candidate BL Lacs. Ten 2BL objects remain 
unclassified; based on contamination levels evident from the 
VLT spectra we may expect one or two of these unclassi- 
fied objects to be a candidate BL Lac, giving a sample of 
~10 optically selected BL Lac objects. Further high signal- 
to- noise, high resolution spectroscopy and/or NIR imaging 
is needed for these ten objects to determine their identity. 

Optical selection of BL Lac candidates can be an ef- 
fective method of identifying candidate objects, providing 
the spectra are flux calibrated and of high enough signal- 
to-noise and resolution to eliminate contamination by weak- 
lined AGN. Overcoming contamination by cool, featureless 
WDs is more difficult (and hence time-consuming), and re- 
lies on accurate proper motion data together with colour 
selection strategies in order to eliminate the bulk of ther- 
mal objects in the initial selection procedure. However some 
overlap in colourspace between WDs and QSOs remains in 
the optical regime. NIR imaging or spectroscopy of radio un- 
detected candidate BL Lacs is thus necessary to distinguish 
between nonthermal objects and cool DC WDs. 

For the 2BL sample the selection criteria of the original 
2QZ and 6QZ catalogues were less than optimal for finding 
candidate BL Lac objects - the bright limit for 2QZ sources 
of 6,j > 18.25, and the requirement that objects be classed 
as point sources, will have eliminated many nearby BL Lac 
objects from the 2QZ catalogue (indeed six BL Lacs were 
independently identified in the 2dF Galaxy Redshift Sur- 
vey (Colless et al. 2001), Lewis, I. J., private communication). 
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Figure 8. Top: a TO — a ox diagram of 2BL BL Lacs; Solid squares 
denote objects with both radio and X-ray fluxes, while triangles 
are those with no X-ray flux. Arrows indicate upper limits of 
fx = 8 X 10~ 14 erg s _1 cm~ 2 . The open squares are the two 
radio-quiet candidate BL Lacs. Bottom: a IO — o x diagram from 
Paper 1 (Fig 14) showing 47 radio-quiet candidate BL Lacs. 



However, from a search of the online Veron & Veron-Cetty 
(2000) catalogue, conducted once the original 2BL sample 
had been selected, we determined that only two possible BL 
Lac objects (denoted BL? by Veron & Veron-Cetty) were 
missed as a result of the 2QZ's initial colour (UV excess) 
selection criteria. 

While there is evidently not a large population of radio- 
quiet BL Lacs, optical selection does probe to lower a ro val- 
ues than do X-ray and radio-selected samples; inclusion of 
these faint candidate BL Lacs is clearly necessary if evolu- 
tionary trends for BL Lac objects are to be correctly inter- 
preted. 
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APPENDIX A: SDSS SPECTROSCOPY OF 2BL 
OBJECTS 

APPENDIX B: OPTICAL AND NIR 
PHOTOMETRY OF 2BL OBJECTS 

APPENDIX C: COMPARISON OF 2DF/6DF 
AND VLT SPECTRA 
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Figure Al. SDSS spectra of 2BL objects 
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Figure Bl. Plots of 2BL objects with SDSS optical magnitudes and NIR magnitudes from Kitt Peak, NOT and Calar Alto 
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Figure CI. 2dF spectra (top) and VLT spectrum of the radio-quiet object J003058.2-275629. 
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Figure C2. Top: 6dF and VLT spectrum of the radio-quiet object J022044.7— 30. Centre: 2dF and VLT spectrum of the radio-detected 
object J031056.9-30. Bottom: 2dF and VLT spectrum of the radio-quiet object J103607.4+01. 
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Figure C3. Top: 2dF and VLT spectrum of the radio-quiet object J120015.3+00. Centre: 2dF and VLT spectrum of the radio-quiet 
object J123437.6-01. Bottom: 6dF spectra and VLT spectrum of the radio-quiet object J220222.7-28. 
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Figure C4. Top: 2dF and VLT spectrum of the radio-quiet object J220850.0— 30. Centre: 2dF and VLT spectrum of the radio-quiet 
object J221105.2-28. Bottom: 2dF and VLT spectrum of the radio-quiet object J221450.1-29. 
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Figure C5. Top: 2dF and VLT spectrum of the radio-quiet object J230306.0— 31. Centre and bottom: 6dF spectra and VLT spectrum 
of the radio-quiet object J233508.1-28 
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